close all;
clear all;
clc;

Nx=10;
Ny=10;
A=zeros(Nx,Ny);
xa = linspace(-1,1,Nx);
ya = linspace(-1,1,Ny);
[X,Y] = meshgrid(xa,ya);
A=(X.^2+Y.^2)<=0.75;



